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ABSTRACT 



Massive binary stars may constitute a substantial fraction of progenitors to 
supernovae and 7-ray bursts, and the distribution of their orbital characteristics 
holds clues to the formation process of massive stars. As a contribution to secur- 
ing statistics on OB-type binaries, we report the discovery and orbital parameters 
for five new systems as part of the Cygnus 0B2 Radial Velocity Survey. Four 
of the new systems (MT070, MT174, MT267, and MT734=VI Cygni #11) are 
single- lined spectroscopic binaries while one (MT103) is a double-lined system 
(B1V+B2V). MT070 is noteworthy as the longest period system yet measured in 
Cyg 0B2, with P=6.2 yr. The other four systems have periods ranging between 
4 and 73 days. MT174 is noteworthy for having a probable mass ratio q < 0.1, 
making it a candidate progenitor to a low-mass X-ray binary. These measure- 
ments bring the total number of massive binaries in Cyg OB2 to 25, the most 
currently known in any single cluster or association. 

Subject headings: Techniques: radial velocities — (Stars:) binaries: general - 
(Stars:) binaries: spectroscopic — (Stars:) binaries: {including multiple) close 
— Stars: early-type — Stars: kinematics and dynamics 



1. Introduction 



Massive stars produce some of nature's most ener getic phenomena as they end their 
lives in spectacular supernova or 7-ray burst explosions (IBurrows et al.lll995l ; IWoosley et al. 
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19931 ; IWoosley fc Bloomll2006l ). While the post-main-sequence evolutionary path of a single 
massive star is still not settled (e.g., the sequence through re d/blue supergiant, Wolf-R ayet, 
Luminous Blue Variable phases for stars of various masses; iMaeder fc Meynetl 120001 ) . the 
possible evolutionary paths of close massive bi nary systems are even less certain. Never- 



theless, close massiv e binaries appear common (IGarmany et al.l Il980l ; ISana fc Evand 12011 



Kiminki et al.l 120121 ). and secondary stars in such systems have been implicated in remov- 
ing hydrogen envelopes from massive primaries prior to core collapse, resulting in a sub- 
stantial fraction (30-40% ) of Type Ib/c super novae and 7-ray bursts (jlzzard et al.l 12004 ; 
Kobulnicky &: Fryer! 120071 ; lEldridge et al.l 120081 ). Close massive binaries may also produce 
the observed population of high-velocity "runaway" stars either through ejection of one 
component after the supernova ex plosion or during multi-body gravitational interactions 
flGies fc Boltor]ll986HBlaauJll96lh . 



Papers I-IV in this series (IKiminki et al.ll2007l iKiminki et al.ll2008[ Kiminki et al.ll2009 



Kiminki et al.ll2012l ) describe prior results of the Cygnus OB2 Radial Velocity Survey in- 
tended to measure the massive binary characteristics (i.e., binary fraction, distribution of 
periods, mass ratios, eccentricities) for a large number (114) massive stars in a single clus- 
ter/association having a common formation environment. In particular, Paper IV (Table 6) 
su mmarizes orbital elements for the 20 massive binaries with measured parameters. Paper 
V ( IKiminki fc Kobulnickyll2012[ ) uses these previously published data and the new data pre- 
sented herein (a total of 25 measured systems) to infer the intrinsic distributions of massive 
binaries, concluding that the fraction of massive stars (defined as B3 and earlier; i.e., su- 
pernova progenitorgj) having companions may be as high as 90%, and that 45% of these 
are likely to interact at some point. Paper V reports an excess of short-period 4-7 day 
systems relative to 7-14 day systems and finds that the distribution of mass ratios is ap- 
proximately flat over the range 0.1< q <1.0. Besides providing the fundamental data for 
modeling the frequencies of energetic phenomena, these types of statistics help place con- 
straints on theoretical frameworks for the formation of massive stars which remain under 



debate flKrumholz et al.ll2010l : ISmith et al.ll20091 ). 



In this sixth paper of the series we present orbital solutions for five additional mas- 
sive spectroscopic binaries in Cyg OB2 — one double-lined (SB2) and four single-lined (SB1) 
systems — using spectroscopic data obtained primar ily during 2010-2011, but ut ilizing some 
data as early as 1999. Using the nomenclature of iMassey fc Thompson! (jl99ll ). these sys- 
tems are MT103 (SB2), MT070, MT174, MT267, and MT734 (SBls). Figure □ displays a 
three-color representation of the Cyg OB2 vicinity, where blue/green/red depict the Palo- 



1 For a typical initial mass function (|Salpetenll955j ). B stars outnumber O stars 3:1 as supernova progen- 
itors. 
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mar Sky Survey R, the Spitzer 4.5 /im, and the Spitzer 8.0 /im bands, respectively. White 
labels denote previously known binary systems, while magenta label s highlight the newly 
discov ered systems reported herein. Numeration follows t he system of IMassey &: Thompson 
( 119911 ), with "S" additi o nally indicating the numeration of lSchultd (119581 ) and "A" or "B" for 



that of IComeron et al.l (120021 ) . Figure [TJ shows the complex nature of this region, including 
ionized gas (diffuse blue emission tracing Ha within the POSS R band), photo-dissociation 
regions (PDR) at the edges of molecular clouds (diffuse red and green tracing broad emission 
features arising from excited polycyclic aromatic hydrocarbons, PAHs), and stars (blue and 
green point sources). 

After first reviewing the Survey strategy and data, we present the new radial velocity 
data and describe the orbital characteristics of these O and early B systems. We conclude 
by summarizing the parameters of the 25 currently known massive binaries in CygOB2. All 
reported velocities are in the heliocentric frame unless explicitly indicated otherwise. 



2. Spectroscopic Observations, Reductions, and Radial Velocity Measurement 

The Cygnus OB2 Radial Velocity Survey is an ongoing o ptical spectroscopic su r vey o f 



~120 massive primary stars (earlier than B3) drawn from the IMassey fc Thompson! (119911 ) 
photometric study of Cyg OB2. Papers I-IV provide details regarding the variety of obser- 
vatories, dates of observation, and data analysis methods employed dating back to the start 
of the survey in 1999. Observing cadences vary from nightly to as few as several observations 
per year for the survey overall. Here, we repeat essential observational details and report 
on new observations covering the 2010 and 2011 observing seasons at the Wyoming Infrared 
Observatory (WIRO) 2.3 m telescope. Table [TJ provides a log of observation dates and de- 
rived radial velocities for the SB2 MT103 and Table [2] lists observing dates and velocities for 
the four SB1 systems. 

In 2011, the data were obtained using the WIRO Longslit spectrograph with an e2V 
2048 2 ccd as the detector. An 1800 1 mm grating in first order yielded a spectral resolution 
of 1.5 A near 5800 A with a 2"xl00" slit. The spectral coverage was 5250-6750 A. Exposure 
times ranged from 1200 s to 3600 s in multiples of 600 s depending on target brightness, cur- 
rent seeing (1.2-3 arcseconds FWHM) and cloud conditions. Reductions followed standard 
longslit techniques, including flat fielding from dome quartz lamp exposures. Copper-argon 
arc lamp exposures were taken after each star exposure to wavelength calibrate the spectra 
to an rms of 0.03 A (1.5 km s _1 at 5800 A). Multiple exposures were combined yielding final 
S/N ratios typically in excess of 100:1 near 5800 A. Final spectra were Doppler corrected 
to the heliocentric frame. Each spectrum was then shifted by a small additional amount in 
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velocity so that the Na I D AA5890,5996 lines were registered with the mean Na I line wave- 
length across the ensemble of observations. This zero-point correction to each observation 
is needed to account for effects of image wander in the dispersion direction when the stellar 
FWHM of the point spread function was appreciably less than the slit width. Because of 
these inevitable slit-placement effects on the resulting wavelength solutions at the level of 
<10 km s _1 , radial velocity standards were not routinely taken. Relative velocity shifts were 
generally less than 7 km s" 1 . 

In 2010 the WIRO-Longslit data were obtained with the 600 1 mm -1 grating in second 
order covering 4000-5900 A at 2.8 A FWHM resolution near 5800 A. Reductions followed the 
same procedures previously described. The lower spectral resolution of these data limited 
their utility and are seldom used, but we note them here for completeness. 

We measured the radial velocity for each spectrum obtained at WIRO using Gaussian fits 
to the He I A5876 line. Systemic velocities are, therefore, based solely on this line, by adopting 
a rest wavl ength of 5875.69 A measured in model stel lar spectra with static atmospheres 



(TLUSTY; lLanz fc Hubenyll2003uHubeny fc Lanzlll995l ). Stars with strong winds, including 



very early O stars (not present in this particular subsample) and evolved stars (MT267, 
MT734) may exhibit He line centers that are blue shifted with respect to this assumed rest 
wavelength. Our fitting cod^l fixes the Gaussian width and depth to be the mean determined 
from all the spectra, after rejecting outliers, and it solves for the best fitting line center. In the 
case of an SB2 like MT103, the code fits a double Gaussian where the widths and depths have 
been fixed independently using observations obtained near quadrature. In the case of the 
long-period system MT070 (6.2 years) we included data from the WIYN telescope & Hydra 
fiber spectrograph covering the range 3800-4500 A. Therefore, we measured radial velocities 
by fits to the He I A4471 line instead of A5876. Our analysis of Keck spectra of stars having 
minimal radial velocity variations suggest no systematic differences between these two lines. 
However, the velocities from the WIYN/Hydra appear systematically 5-15 km s~ x more 
negative than velocities from WIRO for constant-velocity systems having sufficient data 
from both observatories to make this comparison. Observations of radial velocity standard 
stars using WIYN/Hydra showed excellent agreement with published velocities (Paper I). 
Given the spectral stability of the fiber-fed Hydra bench-mounted spectrograph, we consider 
velocities from this instrument to be more reliable. Accordingly, when WIRO and WIYN 
spectra are combined to compute solutions for the same star (MT070), we shift the WIRO 
velocities by —10 km s -1 . We caution that this may introduce systematic shifts in the 
systemic velocity reported. However, the primary goal of the Cygnus OB2 radial velocity 



2 We use the robust curve fitting algorithm MPFIT as implemented in IDL (|Markwardtl l2009). 
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survey is to obtain orbital parameters for massive binaries, a goal that requires good relative 
radial velocities. Absolute space velocities from which systemic velocities may be inferred are 
of lesser importance. Tables [T] and |5] report the adopted radial velocities and uncertainties 
for each star, along with the orbital phase, 0, from the best-fitting solution, and the observed 
minus computed (O— C) residuals. The final column specifies the observatory used to obtain 
the data. 



3. Orbital Solutions 



We analyzed the CLEANed ([Roberts et al.l 119871 ) power spectrum for each object to 



select likely periods and then examined the folded velocity curve for periods corresponding 
to the strongest peaks in the power spectrum. In most cases the strongest peaks yielded 
clean, convincing phased velocity curves. Secondary peaks and possible aliases could be 
eliminated by visual inspection owing to the much larger dispersion in the data at any given 
phase. We used the binary orbital solution package "binary" by D. Gudehus^l with these 
initial period estimates and the radial velocity data to solve for the full suite of orbital 
parameters and associated uncertainties. Table [3] compiles these best-fitting parameters and 
uncertainties for each object. Listed within the table are the period in days (P), eccentricity 
of the orbit (e), longitude of periastron in degrees (a;), systemic radial velocity (7), epoch 
of periastron (T ), primary and secondary (for MT103) projected semi-amplitudes (K 1 & 
K 2 ), spectral classifications from this survey (S.C.i & S.C2), estimates of the mass ratio 
(q), inclination (i), semi-major axis (a), and reduced chi squared values of the best fitting 
solution. The ensuing subsections elaborate on the detailed orbital solutions for the five 
massive systems. 



3.1. MT070 



MT070 has V = 12.99 with a systemic velocity of —6.4 km/s, as reported by lKiminki et al 



(120071 ). Its B — V color of 2.10 mag makes it among the reddest early-type stars in Cyg OB2. 
Figure [T] show that MT070 lies along the eastern edge of the Association (upper right) in 
an apparently unremarkable region mostly devoid of Ha and PAH emission. However, there 
are filaments of 8 fim & 24 /xm emission and 12 CO (C. Brunt, private communication) at 
Cyg OB2 radial velocities cutting across this region, suggesting the likelihood of elevated 
dust extinction that could explain the redder color of MT070 relative to the other Cyg OB2 



3 http://www. chara.gsu.edu/ gudehus/binary.html 
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members. 

Figure |5] displays an average of ten WIRO spectra of MT070 (top) covering the red spec- 
tral range and six WIYN spectra (bottom) covering the blue spectral range. Labels denote 
key spectral features. The composite WIYN spectrum for MT070 shows He II A4200:He I 
A4144, one of the primary temperature-sensitive line ratios in this wavelength r egime, in a 



ratio of 1:1, indicating the temperature class of th e primary is ap p roxim ately Q9 (jSota et al. 



201ll ). in agreement with the initial estimate of iKiminki et al.l ( 120071 ) . Additionally, the 



luminosity-sensitive ratio Si IV A4089:He I A40 26 appears in a r atio of approximately 3:4 



indicating an evolved luminosity class of III — II (ISota et al.ll201ll ). In the composite WIRO 
spectrum, C III A5696 appears to be slightly in emission, O III A5592:He II A5411 appears in 
a ratio of approximately 3:5, and C I V A5801:He II A 5411 appears in a ratio of approximately 



1:2, also indicating an evolved type (jWalbornl 120091 ). We therefore adopt the spectral type 



of 09III for the primary of MT070. This is consistent with the relative faintness (V =12.99) 
but the larger extinction of AAy=+1.8 mag (inferred from B — V=2.10 mag) relative the 
average Cyg OB2 member (B — V~ 1.5 mag). 

Our spectroscopic dataset on MT070 includes 28 epochs covering the period 1999- 
2011 from the Keck (HIRES spectrograph; 2 measurements), WIYN (Hydra spectrograph; 5 
measurements), and WIRO (WIRO-Longslit spectrograph; 21 measurements) observatories. 
Our time-series analysis of the MT070 spectra shows a peak in the power spectrum near 
6 years, along with several secondary peaks at similarly long periods between 5-8 years. 
However, the best-fitting orbital solutions invariably converged to a period of 6.19±0.13 
years, which yielded Xr e d=0.68, considerably smaller than other candidate periods. Figure [3] 
shows the data and best fitting orbital solution. Filled circles, triangles, and squares denote 
the Keck, WIYN, and WIRO data, respectively. This solution specifies a systemtic velocity 
7=— 13±1 km s -1 , eccentricity e=0.34±0.11, and a (projected) primary velocity amplitude 
fr!=9±l km s _1 . Table [3] lists the full orbital parameters for MT070 and other systems. 
The four data points near phase = 0.9 are from WIRO in 2010, and, despite the larger 
uncertainties, these are included to constrain a portion of the orbit otherwise lacking data. 
The radial velocities reported in Table [2] for data from WIRO have been corrected by —10 
km s _1 to bring it into better agreement with the WIYN data. 



If we adopt a primary mass of Mi =22 M based on the 09III spectral type (IMartins et al. 



20051 ). we can begin to place constraints on the secondary mass M2 and on the inclination 
of the system. If the inclination of the system were as low as i=18°, the implied secondary 
mass approaches that of the primary. Making the motivated assumption that the secondary 
is on the main sequence, the luminosity of such a star would be smaller than that of the 
primary, and indeed, we do not see spectral features from the secondary. Therefore, we can 
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limit the inclination to i > 18°, making the mass ratio q = M2/M1 < 1 and semi- major axis 
a < 10.9 AU. Stronger constraints are not possible because of the low velocity amplitude of 
the system relative to the spectral resolution of the data. On the other extreme, if i = 90°, 
the implied secondary mass is 5.0 (approximately a B6V), (7=0.22, and a=9.3 A U. Pho- 
tometric surveys of Cyg OB2 (IPigulski fc Kolaczkowskilll998t iHenderson et al.ll201ll ) do not 
list MT070 as an eclipsing system, but, given the long period, the lack of observed eclipses 
does not place useful constraints on the inclination of the system. 

For completeness, we also report briefly the orbital solution that would have resulted 
without applying the -10 km s _1 systematic velocity correction to the WIRO data. In this 
case, we find a a similar period but a much more eccentric system: P = 5.79 ± 0.11 yr, 
e = 0.71 ± 0.22, a systemic velocity 7 = — 1 ± 4 km s _1 , and a velocity amplitude Ki=24A 
km s _1 . The xl e d * s a ^ so considerably larger at 1.12. Consequently, we adopt the 6.19 yr 
period and solution described previously. In summary, MT070 is a single-lined long-period 
binary with P = 6.19 yr, e ~ 0.34, 0.22 < q < 1, and 18° < i < 90°. 



3.2. MT103 

MT103 lies just 195" from MT070, and it is similarly red, with B - V=2.00. It is the 
only double-lined binary reported in this work, but many SB2 systems within Cyg OB2 
were described in Papers II-IV. The components are significantly blended, except for the 
times very near quadrature. We determine velocities using double Guassian fits to the He I 
A5876 line for the 12 WIRO data (all from 2011) and A4471 for the three WIYN data (two 
from 2001 and one from 2006). We chose to constrain the width of the lines (we used a 
FWHM of 3.0 A, estimated from the most deblended spectra) and the ratio of component 
line depths to be near 2.5:1, within 20%. The power spectrum of the primary exhibited 
a singular strong peak near 22 days with no credible aliases. The power spectrum of the 
secondary was consistent with this value but considerably more noisy owing to larger velocity 
uncertainties. As with MT070, we applied a —10 km s -1 shift to all of the WIRO data to 
improve agreement with the WIYN velocity scale. 

Figure H] shows the best-fitting orbital solution performed jointly for the primary and 
secondary stars. The period of 22.104±0.002 days is quite secure, anchored by the large 
time baseline of the combined WIYN and WIRO data. The eccentricity of e = 0.32 ± 0.05 
is likely to be less secure, owing to the difficulty in resolving the blended spectral features. 
The velocity semi-amplitudes of f^i=71±4 km s _1 and i^2=100±8 km s _1 for the primary 
and secondary, respectively, indicate a mass ratio of 0.71±0.10. The velocity data for the 
secondary has large uncertainties, particularly near <fr =0.6-1 where the deblending technique 
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is only marginally effective. Table |3] includes the full list of orbital elements. 



Figure \5\ (left panel) displays a portion of the 12 spectra around He I A5876 from the 
2011 WIRO observing campaign, ordered by phase. The right panel shows the portion 
of the spectra surrounding Ha. The separation of the components begins near 0=0.16 as 
the primary becomes blueshifted, and the line separation is well-developed at the level of 
~200 km s _1 at 0.38 < < 0.53. Near = 0.72 the primary becomes redshifted and the 
secondary's line is visible on the blueshifted wing until ~ 0.95 or later. 



Kiminki et al.l (120071 ) list MT103 as a B1V. A qualit ative analysis of the com posite of 3 



WIYN spectra covering the blue wavelength range, using iGray fc Corballyj ( 120091 ). supports 
this classification. In the composite spectrum, the Si IV AA4089,4116 lines are absent. This 
suggests that the primary component is a Bl or later. Mg II A4481 is still very weak relative 
to He I A4471, and Si II A4128-30 is also very weak, suggesting that it is earlier than B3. The 
absence of O II A4070, O II A4348, O II A4416 and N II A3995, which we acknowledge may be 
lost in the noise of the spectral continua, indicates that the primary is not evolved, and lies 
in the range B1-2V. We adopt B1V for the primary. The most deblended spectra indicate 
that the secondary component's He lines are similar to the primary's, but slightly weaker. 
Given this and that the calculated mass ratio for this system is 0.71, we adopt B2 as the 
likely spectral type for the secondary. Our attempt at a spectral synthesis of a B1V (14.2 
M Q ; My=— 3.2jj and a B2V (10.9 M Q ; My=— 2.45) provides reasonable agreement with 
the most deblended spectra and is consistent with the measured mass ratio. This suggests 
a luminosity ratio of about 2:1 and a mass ratio consistent with the measured 0.71 value. 
These masses imply an inclination of i ~ 50° in order to match the observed radial velocities. 
The resulting semi-major axis of the systems is then a ~0.45 AU. 

The measured systemic velocity of — 28±3 km s _1 is noteworthy for being about 15 
km s -1 more negativ e than the mean of Vf ie u =—13 km s _1 for massive stars in Cyg OB2 
( IKiminki et al.ll2007l ). Without the 10 km s -1 velocity offset applied to the WIRO data, the 
systemic velocity for MT103 would be in better agreement with the rest of Cyg OB2 massive 
stars. We note here that the uncertainty on the systemic velocity is likely to be somewhat 
larger than the formal 3 km s -1 error quoted and that we do not consider the difference from 
— 13 km s _1 to be significant, given possible systematic velocity uncertainties at the 10-12 
km s _1 level. 



Interpolated from Drilling fc Landolt ( 2000l) . 
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3.3. MT174 



Kiminki et al.l (120071 ) classify MT174 as a B2IV, while Massev fc Thompsonl (Il991f ) es- 
timated B2V. Figure |6] displays an average spectrum of MT174 using ten 2011 WIRO data 
(top) covering the red spectral range and using seven WIYN spectra (bottom) covering 
the blue spectral range. Labels mark key spectral features. The prim ary temperature- 
sensit ive line ratio in this wavelength regime is Si IV A4089:Si III A4552 iGray fc Corbally 
(120091 ). Unfortunately, the composite WIYN spectrum does not cover Si III A4552. How- 
ever, the spectrum does show very weak Si IV A4089 and no Si IV A4116, suggesting a 
temperature class of B1-B2. The He I A4471:Mg II A4481 temperature-sensitive ratio, 
which becomes useful around B2-B3 is ~1:5 and indicates a temperature class of B2 or 



earlier (IWalborn fc Fitzpatrick! Il990f ). Additionally, He I A4009 and He I A4026 are quite 
stro ng and in a ratio of 1:2 respec tively, also indicating a temperature class of approximately 
B2 ( IWalborn fc Fitzpatrick! Il990l ). The composite WIYN and WIRO spectra show narrow 
Balmer line widths. In conjunction with the weak, but clearly present N II A3995 and O II 
AA4070,4076, this suggests a luminosity cl a ss of H I. Finally, a direct comparison to HD35468 
(B2III) from the IWalborn fc Fitzpatrick! (jl990l ) digital atlas provides the best agreement 
with both the composite WIYN and WIRO spectra. Therefore, we adopt a spectral type of 
B2III for the primary component of MT174. At V=12.55 and B-V=1.21, it is not unusually 
colored for early type stars in Cyg OB2, but it is 1-2 magnitudes brighter than other B2 
stars in our CygOB2 sample, consistent with an evolved luminosity class. 

Our 10 WIRO spectra from 2011 reveal that this is a single-lined spectroscopic binary 
with a period of 4.536±0.020 d, e = 0.53 ± 0.13, and Ki=9±l km s -1 . Figure [7] shows the 
data and best fitting orbital solution, which has Xred = ®-7^- The systemic velocity of +11 ±2 
km s _1 is considerably more p ositive than the —13 km s _1 average of probable Cyg OB2 
members (IKiminki et al.l 120071 ) . This may be the result of systematic differences between 
the WIRO data and that from other observatories. Our attempts to include the seven 
measurements from WIYN over the period 2001 - 2006 resulted in poor or non- convergent 

solutions, even after applying a —10 20 km s" 1 offset to the WIRO measurements. This 

may indicate the presence of additional kinematic variations in the MT174 system, possibly 
the result of the gravitational influence of a third body on timescales of many months or years. 
Unfortunately, our most reliable data do not have sufficient time coverage to constrain the 
nature of possible additional velocity perturbations in MT174. Given that the mean velocity 
measured in our WIYN spectra is Vh e i io =—Q km s _1 compared to +12 for the WIRO data, 
we apply a —18 km s _1 zero-point shift to our final orbital solution, summarized in Table [3] 
and the radial velocities in Table [21 The resulting systemic velocity of -6 km s _1 is more 
consistent with the ensemble of Cyg OB2 stars. 
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If we adopt a primary mass appropriate to a B2III (M ~10.9 M ; i?=10 Rq), the implied 
secondary mass is M 2 =0.35 M for i =90° and M 2 =1.0 M for i =20°. For inclinations as low 
as 3° M 2 approaches M 1; a lowe r bound imposed by th e lack of secondary spectral features. 
The absence of known eclipses ([Henderson et al.ll201ll ). coupled with the narrow linewidth 
of ~1 A FWHM for He I A4471 in our single 1999 Keck spectrum, is consistent with a small 
projected rotational velocity, which would imply a small inclination angle if the components' 
spin axes are aligned with the orbital axis. (The small linewidth is also consistent with the 
evolved luminosity classification.) The inclination is, effectively, unconstrained, and may lie 
in the range 3° < i < 87°. For all inclinations greater than about 15 degrees, the secondary 
is M<1.4 M , making MT174 one of the most extreme mass ratio (q < 0.13) systems yet 
uncovered in Cyg OB2. 

MT174 is an interesting system in that it is the first binary in the Cyg OB2 radial 
velocity survey to contain a massive star and a quasi solar-mass star. As such, it is a 
potential progenitor of a low- mass X-ray binary system (LMXRB). The distribution of mass 
ratios among Cyg OB2 stars is approximat ely flat, and some 45% of mas sive binary systems 
have components capable of mass transfer, (jKiminki fc Kobulnickyll20f 21 ). meaning that ~5- 
10% of massive binaries may be progenitors of low-mass X-ray binaries. However, the relative 
paucity of observed LMXRBs in nature implies that the formation channels for such objects 
ar e narrow and require sp ecial conditions during post-main-sequence evolution, as discussed 
by Ivan den Heuvell ( 119831 ) . 



3.4. MT267 (All) 

Figure M shows the 25 WIRO data, all from 2011 June - 201 1 November, and th e 



orbital solution for the SB1 system MT267 (All in the notation of IComeron et al.ll2002l ). 
The period of P =15.511±0.056 d, eccentricity of e =0.21±0.07, and velocity amplitude 
A' 1 =24±2 km s _1 are well constrained by the data. However, the xled °f 1-59 is larger than 
for most other systems, suggesting the possibility of photospheric line variations (common for 
evolved massive stars) or an additional dynamical influence in the system. Our examination 
of the observed minus computed (O-C) values in Table [2] reveals a weak correlation between 
Julian date and O-C, but additional data will be needed to establish the nature of additional 
long-term velocity variations, such as those expected from a third body. The resulting 
systemic velocity of 7=— 13±1 km s -1 is consistent with other CygOB2 stars, and, given 
the lack of data from other observatories, we do not apply any systematic correction to the 
velocities plotted in Figure |8] or the data in Table |2j The full suite of orbital parameters 
appears in Table [31 
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Figure |9] is an average spectrum of MT267 covering the wavelength range 5400-6700 A 
from our WIRO data. Labels mark key spectral atmospheric and interstellar medium fea- 
tures. The composite spectrum exhibi ts a He II A54 11:He I A5876 ratio of near unity, 
indicating a temperature class near 07 (jWalbornl Il980l ) . O III A5592 is present at r oughly 
1:3 with He I A5876 and also agrees with a late-0 temperature class (jWalbornl Il980l ). C III 
A5696 is prominent in emission. C IV AA5801-5812 and O III A5592 are present at moderate 
strength and relatively br oad, and Ha appears to be partially filled in by emission. This all 
indicates an evolved star (jWalbornl 120091 Il980l ). Ha shows variable emission, which is also 
indicative of an evolved star or an interacting binary. Figure [10] displays a sequence of Ha 
(right) and He I A5876 (left) spectra, ordered by phase. Changes in the Ha line profile are 
evident, but He I appears u nchanged, except for the periodic velocity variations. Comparing 
spectral standards from the IJacoby et al.l (119841 ) digital atlas with a composite spectrum of 
this star, we find that an 07-08 best fits the temperature of primary component, with the 
luminosity class in the range III — I. This also agrees with our photometric assessm ent of this 
star and is consistent with the spectral type estimated by iNegueruela et al.l (120081 ) who give 
07.5Ib-II(f). 



If we adopt a mass of 32 M appropriate for an 07.5111 - I (jMartins et al.l 120051 ) then 
the implied secondary mass ranges between 3.0 M for an inclination of i = 90° and 30 
M Q for i = 8° if the secondary approaches the mass of the primary. Given the evolved, 
luminous nature of the primary, which we assume dominates the luminosity of the system, 
is it conceivable that the secondary may be as massi ye as the primary but not visible in 
the integrated spectrum. The lack of known eclipses (jHenderson et al.l 1201 ll ) provides only 
a weak constraint that the inclination is smaller than 81°. Thus, we conclude that the mass 
ratio of the system is 0.09 < q < 1 with semi major a xis 0.48 > a > 0.39 AU. The presence of 
X-rays and X-ray variability in MT267 (jRauwll201ll ) is consistent with a wind-wind collision 
in a close, interacting binary. If the wind collision scenario is correct for the origin of the 
X-rays, then the secondary component is probably toward the massive end of our allowed 
range, implying i < 13° in order to make the mass of the secondary equivalent to a B0V or 
earlier. 

MT267 has an unusual ly red B — V color (2 . 19 m ap;) compared to other massive stars 
in Cyg OB2 (<B-V>~1.4; iMassey &: Thompson! 1 19911 ). meaning a higher foreground dust 
column density toward this star. Figure ITT1 shows an infrared view of the field around MT267. 
The left panel is a three-color image with blue/green/red representing the POSS R-band, the 
Spitzer 4.5 /zm band, and Spitzer 24 /im band. In the right panel blue/green/red represent 
the UKIDSS JHK images. Labels mark the locations of MT267, and the well-known, heavily 
reddened mid-B supergiant Cyg OB2 #12 (MT304). Interestingly, these two evolved massive 
stars lie within a radius of 53" from each other. Clearly evident in the left panel is the infrared 
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dar k cloud adjac e nt to both objects. This is the Bok globule noted by lScappini et al.l (120071 ) 
and iPoppi et al.l ( 120101 ) . Green contours show the 12 CO line intensity integrated over the 
velocity range Vlsr=9.7 - 11.7 km s^ 1 , corresponding to V Q =23.7 - 25.7 km s _1 |§ Contour 
levels show 13, 25, 38, and 50 K km s -1 intensities corresponding to column densities of 
N(H2)=3.9, 7.5, 11.4, and 15.0 xlO 21 cm" 2 , using the mean relation bet ween molecular 
colum n density and 12 CO integrated i ntensity a = 3 .0 x 1 20 cm 2 / (K km s _1 ) ( jSolomon et al. 



19871 ). Given the conversion factor of iBohlin et a. dl978)B the contours correspond roughly 
to extinctions A v =0.67, 1.3, 2.0, 2.6 mag, assuming R v = 3.1. These are probably lower 
limits given that 12 CO becomes optically thick and/or freezes onto dust grains in the cold 
cores of molecular clouds. If this Bok globule/molecular cloud lies in the foreground to 
MT267 (and Cyg OB2#12), it is conceivably responsible for the higher extinction recorded 
to these objects. However, the very large extinction of Cyg OB2#12 (Ay ~12 mag) cannot 
be solely attributed to this cloud. Given the large radial velocity difference between this 
cloud and Cyg OB2, it is probable that this cloud lies substantially in the foreground. The 
lack of PAH emission around the borders of this cloud further supports the conjecture that 
it lies at some distance from the massive ionizing stars. 



3.5. MT734 (Schulte #11) 

MT734 (VI Cyg #11 /Schult e 11) is a well-known emission-line star of type 051(f), as 
classified by iHerrero et al.l (120021 ) . Figure [1] illustrates that MT734 lies along the western 
boundary of the Association in an area apparently devoid of Ha and PAH emission, although 
it is surrounded by molecular clouds. Our 39 spectra, taken with the WIRO Longlit spectro- 
graph over the period 2007-2011, show that this system is a single-lined spectroscopic binary. 
We find a period of 72.43±0.07 d, e =0.50±0.06, and A^=26±l km s -1 . Figure Q2] depicts 
our best orbital solution for this system. MT734 possesses a systemic velocity of — 31±1 
km s -1 which is more negative than most members of the Association. We have not applied 
any systematic velocity corrections to the WIRO data. The comparatively blue-shifted line 
centers may be a result of line formation in the upper atmosphere/outflowing wind in this 
very luminous supergiant. 



If we adopt a mass of 51 M Q appropriate to an 051 (R = 19 R & ) from iMartins et al. 

( 120051 ). then an inclination of 90° yields a secondary mass of M 2 =7.7 M , g=0.15, and 



5 The conversion from the local standard of rest (LSR) to heliocentric reference frame is Vq 
km s _1 in the direction of CygOB2. 



Vlsr +14 



6 = < E(B - V) > (mag) 



N(H + H 2 ) cm 2 / 5.8 x 10 21 (cm 2 mag" 1 
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a=1.31 AU. Assuming i=13° implies that the seondary mass approaches that of the primary, 
q ~ 1 and a=1.51 AU. The lack of known eclipses does not place helpful constraints on the 
i nclination. MT734 is listed as a 1400 MHz radio continuum source with S , i 4 oo=0.5±0.3 mJy 
( ISetia Gunawan et al.ll2003l ). This could result from the coliding winds of two massive stars. 

Figure [TBI shows the spectral time sequence for MT734, ordered by phase, in the spectral 
region around He I A5876 (left panel) and Ha (right panel). The He I line shows a smooth 
progression of central wavelength with orbital phase, while the Ha line is broad, exhibits 
emission, and has irregular line profiles, typical of supergiants. 



4. Discussion and Summary 

We have presented orbital solutions for five additional massive binaries that are probable 
members of the Cygnus OB2 Association, bringing the total number of multiple systems 
therein to 25. The vast majority of these have complete orbital solutions, making this the 
largest collection of m a ssive binary statistics in any single cluster or association to date. 



Kiminki &: Kobulnickyi ( 120121 ) have used these data, along with previously known systems 



documented in Papers I-IV, to infer the underlying distrib ution of periods, mass ratios 



and eccentricities for Cyg OB2 as a whole. The analysis of iKiminki &: Kobulnickyi (120121 ) 
suggests that the known list of 25 binaries out of 114 Cyg OB2 systems surveyed is reasonably 
complete at periods less than about 30 days, but that we should expect to discover an 
additional 10-15 binaries having longer periods as the Cygnus OB2 Radial Velocity Survey 
continues. 

Noteworthy among the new list of massive binaries is the discovery of a probable progen- 
itor of a low-mass X-ray binary system (MT174). Also noteworthy is the long-period system 
MT070 with a period of 6.2 years. Both of these low-amplitude systems (Ki~9 km s _1 ) 
attest to the capabilities of the Survey to achieve orbital solutions for either low-inclination 
or low-g systems. 



We are grateful to Gregor Rauw who made us aware of the possible binary nature of 
MT267. We thank the time allocation committees of the Lick, Keck, WIYN, and WIRO 
observatories for granting us observing time and making this project possible. The efforts 
of WIRO staff James Weger and Jerry Bucher made this science possible. We acknowledge 
continued support from the National Science Foundation through Research Experience for 
Undergraduates (REU) program grant AST 03-53760, through grant AST 03-07778, and 
through grant AST 09-08239, and the support of the Wyoming NASA Space Grant Consor- 
tium through grant #NNX10A095H. 
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Table 1. Ephemeris for MT103 







v rl 




V r2 


2 -C* 2 




Date (HJD-2,400,000) 




(km s- 1 ) 


(km s- 1 ) 


(km 




(km s- 1 ) 


Observatory 


MT103 


52146.75 


0.660 


-12.0 ( 5.0 ) 


1.5 


-179.0 ( 


50.0 ) 


-128.6 


WIYN 


52162.75 


0.383 


-102.9 ( 8.0 ) 


-6.5 


28.0 ( 


19.0 ) 


-37.4 


WIYN 


53989.75 


0.036 


18.0 ( 9.0 ) 


7.4 


-186.0 ( 


33.0 ) 


-102.0 


WIYN 


54696.75 


0.031 


7.0 ( 6.0 ) 


-4.3 


-117.0 ( 


23.0 ) 


-32.0 


WIRO 


54697.75 


0.075 


6.0 ( 4.0 ) 


1.2 


-142.0 ( 


45.0 ) 


-66.1 


WIRO 


55779.75 


0.018 


0.0 ( 9.0 ) 


-12.8 


-99.0 ( 


50.0 ) 


-11.9 


WIRO 


55791.50 


0.551 


-55.0 ( 7.4 ) 


27.5 


41.0 ( 


13.0 ) 


-5.0 


WIRO 


55805.50 


0.186 


-16.8 ( 12.9 ) 


2.2 


-55.0 ( 


14.0 ) 


-12.2 


WIRO 


55832.50 


0.403 


-107.0 ( 4.4 ) 


-1.7 


71.1 


( 9-0 ) 


-6.8 


WIRO 


55834.50 


0.495 


-117.7 ( 6.7 ) 


-2.0 


80.5 ( 


12.0 ) 


-11.8 


WIRO 


55855.25 


0.441 


-113.4 ( 6.9 ) 


5.0 


90.0 ( 


11.0 ) 


-6.1 


WIRO 


55857.50 


0.535 


-99.0 ( 5.1 ) 


-4.4 


84.7 ( 


11.0 ) 


21.8 


WIRO 


55866.50 


0.941 


32.5 ( 7.6 ) 


12.6 


-76.0 ( 


34.0 ) 


21.0 


WIRO 


55901.00 


0.511 


-107.0 ( 5.0 ) 


2.1 


101.0 ( 


13.0 ) 


17.9 


WIRO 


55906.25 


0.744 


16.0 ( 11.5 ) 


4.3 


-92.8 ( 


38.5 ) 




WIRO 
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Table 2. Ephemerides for MT070, MT174, MT267, & MT734 





V rl 


O1-C1 


Date (HJD-2,400,000) < 


t> (km s- 1 ) 


(km s _1 ) Observatory 


MT070 



Ol4D ( .UU 


U.1O0 


1(1 7 

-iy. / 


) 


-1.0 




OloUO.UU 


U.olO 


-O. { 


l.O ) 


1 
A. 1 




oz14d. (O 


0.465 


-7.3 


[ 5-4 ) 


0. / 


Wl Y IN 


col fi1 vc; 


n A TO 

U.4fz 


Q Q 

-o.o 


1 A A \ 

4.4 ) 


-u.y 


Wl I IN 


OZWZ. ( O 


U.4fz 


inn 
-1U.U 


4.0 ) 


O 1 
-Z. 1 


Wl I IN 


c^QQ/in trn 
Doo4U.0U 


u.yy4 


Oft A 


O.O J 


O A 


Wl i IN 




281 


-18.0 


9.9 ) 


-5.5 


WIYN 


54399.75 


0.462 


-9.1 


( 4 -8 ) 


-1.1 


WIRO 


54401.75 


0.463 


-16.1 


( 3.3 ) 


-8.1 


WIRO 


54402.75 


0.464 


-6.8 


( 3.6 ) 


1.2 


WIRO 


54410.75 


0.467 


-8.1 


( 4 -2 ) 


-0.2 


WIRO 


54696.75 


0.594 


-5.2 


( 2.9 ) 


1.7 


WIRO 


54697.75 


0.594 


-7.7 


( 3.3 ) 


-0.8 


WIRO 


54698.75 


0.595 


-6.5 


( 2.7 ) 


0.4 


WIRO 


54757.75 


0.621 


-6.3 


( 2.8 ) 


0.6 


WIRO 


55404.75 


0.907 


-9.1 


( 8-6 ) 


7.4 


WIRO 


55422.00 


0.915 


-10.8 


( 9-3 ) 


6.3 


WIRO 


55468.75 


0.935 


-26.1 ( 


12.2 ) 


-7.0 


WIRO 


55491.50 


0.945 


-22.0 


( 3.4 ) 


-1.9 


WIRO 


55716.75 


0.045 


-24.4 


( 2.6 ) 


0.5 


WIRO 


55718.75 


0.046 


-24.1 


( 2.4 ) 


0.8 


WIRO 


55766.75 


0.067 


-23.9 


( 2.2 ) 


0.2 


WIRO 


55805.50 


0.084 


-23.2 


( 2.2 ) 


0.0 


WIRO 


55832.50 


0.096 


-26.1 


( 2.5 ) 


-3.6 


WIRO 


55834.50 


0.097 


-22.2 


( 2.1 ) 


0.3 


WIRO 


55847.50 


0.103 


-22.0 


( 2.1 ) 


0.1 


WIRO 


55855.25 


0.106 


-22.7 


( 2.6 ) 


-0.9 


WIRO 


55866.50 


0.111 


-19.7 


( 2.2 ) 


1.9 


WIRO 



MT174 



55727.00 0.974 7.3 ( 4.3 ) -1.2 WIRO 

55759.00 0.029 16.7 ( 1.9 ) -2.1 WIRO 

55781.00 0.879 6.0 ( 2.9 ) 3.0 WIRO 

55791.50 0.194 19.7 ( 3.1 ) 2.0 WIRO 

55832.50 0.233 15.0 ( 2.6 ) -1.5 WIRO 

55836.50 0.115 20.9 ( 2.5 ) 0.7 WIRO 

55855.50 0.304 13.4 ( 2.8 ) -1.2 WIRO 

55857.50 0.745 4.8 ( 2.6 ) -0.6 WIRO 

55876.25 0.879 1.0 ( 2.6 ) -2.1 WIRO 

55906.25 0.493 11.2 ( 2.5 ) 0.8 WIRO 
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Table 2 — Continued 







V r l 


Oi-Ci 




Date (HJD-2,400,000) 


<$> 


(km 




(km s- 1 ) 


Observatory 


MT267 


^71 ^ 7^ 


0.950 


19.2 


4.2 ) 


4.9 


WIRO 


^71 A 7^ 


0.009 


7.7 


4.0 ) 


-0.9 


WIRO 


^71 ^ 7^ 


0.075 


-7.1 


4.2 ) 


-3.2 


WIRO 


^71 (\ 7^ 


0.138 


-16.2 


4.3 ) 


-0.1 


WIRO 


55717.75 


0.202 


-19.1 


( 4 -3 ) 


6.2 


WIRO 


55718.75 


0.269 


-25.2 


( 4 -3 ) 


5.5 


WIRO 


55726.75 


0.783 


6.1 


( 4.0 ) 


5.7 


WIRO 


55727.75 


0.850 


6.1 


( 3.2 ) 


-2.8 


WIRO 


55737.75 


0.493 


-27.5 


( 4.2 ) 


1.3 


WIRO 


55738.75 


0.558 


-19.3 


( 3.3 ) 


5.1 


WIRO 


55739.50 


0.601 


-24.3 ( 


16.7 ) 


-3.5 


WIRO 


55740.50 


0.665 


-12.7 


( 3.7 ) 


1.6 


WIRO 


55741.50 


0.735 


-11.1 


( 3.3 ) 


-5.1 


WIRO 


55758.75 


0.846 


13.3 


( 3.0 ) 


5.0 


WIRO 


55765.75 


0.294 


-35.8 


( 4.2 ) 


-4.0 


WIRO 


55767.50 


0.406 


-41.0 


( 5.0 ) 


-8.7 


WIRO 


55768.50 


0.472 


-30.6 


( 4.8 ) 


-0.8 


WIRO 


55779.75 


0.192 


-34.7 


( 5.6 ) 


-10.6 


WIRO 


55779.75 


0.200 


-18.3 


( 4.8 ) 


6.8 


WIRO 


55780.75 


0.257 


-35.8 


( 5.9 ) 


-5.8 


WIRO 


55781.75 


0.321 


-37.9 


( 9.1 ) 


-5.4 


WIRO 


55790.50 


0.889 


6.5 


( 3.8 ) 


-6.1 


WIRO 


55834.50 


0.725 


-10.1 


( 4.0 ) 


-2.8 


WIRO 


55857.25 


0.199 


-26.5 


( 4.3 ) 


-1.6 


WIRO 


55866.25 


0.781 


1.5 


( 4.5 ) 


1.5 


WIRO 


MT734 


54347.75 


0.025 


-22.8 ( 


10.0 ) 


-1.1 


WIRO 


54348.75 


0.039 


-17.4 ( 


10.0 ) 


-1.4 


WIRO 


54399.75 


0.743 


-61.2 ( 


16.7 ) 


-16.3 


WIRO 


54402.75 


0.785 


-42.2 


( 8.3 ) 


6.1 


WIRO 


54403.50 


0.795 


-54.9 ( 


10.0 ) 


-5.8 


WIRO 


54403.75 


0.798 


-46.6 ( 


10.0 ) 


2.8 


WIRO 


54405.75 


0.826 


-56.3 ( 


11.1 ) 


-4.5 


WIRO 


54642.75 


0.098 


-8.8 


( 5.3 ) 


-2.6 


WIRO 


54644.00 


0.116 


-6.1 


( 5.6 ) 


-0.1 


WIRO 


54644.75 


0.126 


-6.5 


( 5.3 ) 


-0.4 


WIRO 


54645.75 


0.140 


-5.4 


( 5.6 ) 


1.1 


WIRO 


54646.75 


0.154 


-14.3 


( 5.9 ) 


-7.3 


WIRO 


54647.75 


0.167 


-0.5 


( 5.3 ) 


7.3 


WIRO 


54648.75 


0.181 


-14.7 


( 5.3 ) 


-6.2 


WIRO 


54670.00 


0.475 


-30.0 


( 8.3 ) 


-3.3 


WIRO 
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Table 2 — Continued 



Date (HJD-2,400,000) 




v rl 

(km s _1 ) 


Oi-Ci 
(km s- 1 ) 


Observatory 


54671.75 


0.499 


-33.6 ( 7.7 ) 


-5.3 


WIRO 


54696.75 


0.844 


-62.7 ( 9.1 ) 


-9.4 


WIRO 


54747.75 


0.548 


-33.5 ( 8.3 ) 


-2.2 


WIRO 


54748.75 


0.562 


-26.4 ( 8.3 ) 


5.8 


WIRO 


55404.75 


0.619 


-34.6 ( 12.5 ) 


1.3 


WIRO 


55409.75 


0.688 


-62.8 ( 14.3 ) 


-22.1 


WIRO 


55421.75 


0.854 


-46.7 ( 14.3 ) 


7.4 


WIRO 


55422.75 


0.868 


-50.1 ( 14.3 ) 


5.1 


WIRO 


55423.75 


0.882 


-46.1 ( 14.3 ) 


10.0 


WIRO 


55461.75 


0.406 


-22.4 ( 20.0 ) 


0.2 


WIRO 


55466.75 


0.475 


-24.8 ( 20.0 ) 


2.0 


WIRO 


55498.75 


0.917 


-44.3 ( 20.0 ) 


13.1 


WIRO 


55715.00 


0.903 


-57.9 ( 5.0 ) 


-0.7 


WIRO 


55719.00 


0.958 


-62.8 ( 7.1 ) 


-9.9 


WIRO 


55737.75 


0.217 


-10.2 ( 6.2 ) 


0.5 


WIRO 


55740.75 


0.259 


-7.5 ( 5.9 ) 


5.8 


WIRO 


55780.00 


0.800 


-54.4 ( 6.7 ) 


-4.8 


WIRO 


55790.75 


0.949 


-55.2 ( 5.3 ) 


-0.3 


WIRO 


55794.50 


0.001 


-20.7 ( 12.5 ) 


14.1 


WIRO 


55805.75 


0.156 


-9.2 ( 5.9 ) 


-2.0 


WIRO 


55834.50 


0.553 


-32.5 ( 8.3 ) 


-0.9 


WIRO 


55847.50 


0.732 


-44.0 ( 8.3 ) 


0.1 


WIRO 


55856.25 


0.853 


-51.9 ( 8.3 ) 


2.2 


WIRO 


55857.50 


0.871 


-49.9 ( 7.1 ) 


5.5 


WIRO 
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Table 3. Orbital Elements 



Element 


MT070 


MT103 


MT174 


MT267 


MT734 


P (Days) 


2259±46 


22.104±0.002 


4.536±0.020 


15.511±0.056 


72.43±0.07 


e 


0.34±0.11 


0.32±0.05 


0.53±0.13 


0.21±0.07 


0.50±0.06 


w (dcg) 


157.9±26 


26±12 


282.0±20 


37±19 


262. 6±8 


7 (km s _1 ) 


-13±1 


-28±3 


-6±2 


-13±1 


-31 ±1 


To (HJD-2,400,000) 


53355±134 


55757.2±0.5 


55749.8±0.1 


55745.7±0.8 


55722±1 


K\ (km s- 1 ) 


9±1 


71±4 


9±1 


24±2 


26±1 


K 2 (km s- 1 ) 




100±8 








S. Ci 


09III 


B1V 


B2III 


07.5III-I 


051(f) 


S. C.2 


B: 


B2V 




B - O 


early B - O 


<? 


0.22-1 


0.71±0.10 


0.03-0.9 


0.09-1 


0.15-1 


j (degrees) 


90-18 


~50 


87-3 


81-8 


90-13 


a (AU) 


9.3-10.9 


0.45±0.02 


0.12-0.15 


0.39-0.48 


1.3-1.5 




0.68 


2.41 


0.71 


1.59 


0.54 
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Table 4. OB Binaries in Cyg OB2 



Star 


Type 


S.C. 


(days) 


c 


q 


Ref. 


MT059 


SB1 


08V & B 


4.8527±0.0002 


0.11±0.04 


0.22-0.67 


1 


MT070 


SB1 


09III & B 


6.19 yr 


0.34±0.11 


0.25-1 





MT103 


SB2 


B1V + B2? 


22.104±0.002 


0.32±0.05 


0.71±0.10 





MT145 


SB1 


09III & mid B 


25.140±0.008 


0.291±0.009 


0.26-0.63 


2 


MT174 


SB1 


B2III & ?? 


4.536±0.020 


0.53±0.13 


0.03-0.9 





MT252 


SB2 


B2III & B1V 


18-19 




0.8±0.2 


1 


MT258 


SB1 


08V & B 


14.660±0.002 


0.03±0.05 


0.18-0.89 


1 


MT267 


SB1 


07.5III-I fc O/B 


15.511±0.056 


0.21±0.07 


0.09-1 





MT311 


SB2 


B2V & B3V 


5.7521±0.0002 


0.02±0.01 


0.8±0.1 


3 


MT372 


SB1/EA: 


BOV & B2:V 


2.228 (fixed) 


0.0 (fixed) 


~0.6 


2,4 


MT421 


SB1:/EA 


09V & B9V-A0V 


4.161 




~0. 16-0. 19 


5 


MT429 


SB1/EA 


BOV & GOBI: 


2.9788 (fixed) 


0.38±0.08 


0.057: 


3.5 


MT605 


SB2 


B1V & Bl: 


~4-5 




0.9±0.1 


3 


MT696 


SB2/EW/KE 


09.5V & BOV 


1.4694±0.0002 


0.081±0.005 


0.94±0.03 


3,6 


MT720 


SB2 


B0-B1 & B1-B2 


4.0677±0.0003 


0.35±0.02 


0.80±0.08 


1,3 


MT734 


SB1 


051 & O/early B 


72.43±0.07 


0.50±0.06 


0.15 -1 





MT771 


SB2 


07V & 09V 


2.82105±0.00003 


0.547±0.004 


0.95±0.09 


1.3 


Schultc 3 


SB2/EA: 


06IV: & 09III 


4.7464±0.0002 


0.070±0.009 


0.44±0.08 


1,7 


Schultc 5 


SB2/EB 


07Ianfp & Ofpe/WN9 


6.6 (fixed) 


0.0 (fixed) 


0.28±0.02 


8,9,10,11, 






(& BOV:) 








12,13,14 


Schulte 8a 


SB2 


05.51 & 06: 


21.908 (fixed) 


0.24±0.04 


0.86±0.04 


15,16 


Schulte 9 


SB2 


05: & 06-7: 


2.355 yr 


0.708±0.027 


0.9±0.1 


17,18 


Schulte 73 


SB2 


08III & OSIII 


17.28±0.03 


0.169±0.009 


0.99±0.02 


2 


A36 


SB2/EA 


BOIb & BOIII 


4.674±0.004 


0.10±0.01 


0.70±0.06 


2,19,20 


A45 


SB2 


B0.5V & B2V:-B3V: 


2.884±0.001 


0.273±0.002 


0.46±0.02 


2,20 


B17 


SB2/EB: 


07: & 09: 


4.0217±0.0004 


(fixed) 


0.75 (fixed) 


21,22 



Note. — Photometric types EW/KE, EA, and EB stand for eclipsing system of the W UMa type (ellipsoidal; 
P < 1 day), Algol type (near spherical), and Lyr ty pe (ellipsoidal: P > 1 day) respectively. The mass 
ratio for MT421 is calcu lated using the O star masses of Martins et al. (2005) and interpolated AB masses of 
iDrilling fc Landol j feOOfJ) . 



References. — (0) T h is wo rk 
Pigulski fc Kolaczkowskl Jl998l): (6) iRios 



( 1) Paper II; (2) Paper III ; (3) P aper IV; (4) IWozniak et al.| J2004T) : (5) 
J) | Rios fc peGioia-Eastwoodl J2004); (7) Kiminki fc Kobulnickvl <20ia): (8 ) 
Wilsonl < 19481); f9 l lWiIson fc Abj jl95lh: ( 10'l lMiczaikal <1953l); ( l mWalbornl 19731); (12 ) Contreras et al.l Jl997h: 
(13) | Rauw et al.1 dl99Sh : (14) lHalll <1974h; ( 15) iRomand lll969h; ( 16) | De Becker et al.l |2004h; (17 ) iNaze et alJ 
200sh; (18) iNaze et alj <20ld) ; (19) lOterol i2008al) ; (20) iBansonl j2003 ); (21) IStroud et al.l <20ld) : (22) lOterol 
2008b) 



Fig. 1. — Three-color image of the Cygnus 0B2 region with the Palomar Sky Survey R image 
in blue, the Spitzer 4.5 /zm in green, and Spitzer 8.0 /im in red. White labels denote pre- 
viously known binary systems while magenta labels highlight the newly discovered s ystem s 
reported herein. Generally, numeration follows the system of iMassey fc Thompson! (1l99ll ). 



wi th "S" addition a lly in dicating the numeration of ISchulte rll958h and "A" or "B" for that 
of IComeron et al.l ( 120021 ). The early-B SB2 system A45 lies just off the field of view to the 
upper right of B17. 
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Fig. 2. — Average spectrum of MT070 using data from WIRO (top) and WIYN (bottom). 
Tickmarks and labels mark key spectral features. Mos t of the unlabeled absorpti on features 
in the upper panel are from diffuse interstellar bands (Uenniskens fc Desertlll994l ). 
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Fig. 3. — Radial velocity data and best-fit orbital solution curve for the 6.19 year SB1 
system MT070 (09III:) using 21 measurements from WIRO (filled squares), five from WIYN 
(triangles), and two from Keck (circles). The four points with the largest error bars are from 
the 2010 WIRO campaign which suffered from low spectral resolution. 




Fig. 4. — Radial velocity data and best-fit orbital solution curve for the 22.10 day SB2 
system MT103. Filled symbols denote the primary and open symbols show the secondary. 
Triangles show WIYN data and squares show WIRO data. 
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Fig. 5. — A sequence of the 12 WIRO spectra for MT103, ordered by phase, in the spectral 
region around He I A5876 (left) and Ha (right). The Na D doublet lies to the right of the 
He A5876 line. 
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Fig. 6.— Average spectrum of MT174 (B2III) using data from WIRO (top) and WIYN 
(bottom). 
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Fig. 7.— Orbital solution for the SB1 system MT174 (B2III) using 10 WIRO data from 
2011. 
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Fig. 8.— Orbital solution for the SB1 system MT267 (07.5111 I). 



-32 - 



i i i i | i i i | i i i p 
b-He II III C III C IV He I 

5411 5592 5696 5801,12 5876 



0.7 



_L 



_L 



_L 



Ha 

6563 



He I 

6678 




MT267 



_L 



5400 5600 5800 6000 6200 

Wavelength (Angstroms) 



6400 



6600 



Fig. 9. — Average of 25 WIRO spectra for MT267. Labels mark key spectroscopic features of 
hot stars along with some prominent insters tellar features. Most of the unlabeled absorption 
features are from diffuse interstellar bands (Uenniskens fc Desertlll994l ). 
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Fig. 10. — A sequence of spectra covering Ha and He I A5876 for MT267, ordered by phase. 
This 07III: star exhibits variability in the Ha line profile, but not in the He line profiles. 



Fig. 11. — The field around MT267. The left panel shows the optical/mid-IR view with 
POSS R-band in blue and Spitzer GLIMPSE 8 /im and 24 fim in green and red, respectively. 
Green contours depict the 12 CO globule, which lies coincident with an IR dark cloud. The 
right panel shows the UKIDSS JHK appearance (blue/green/red) of the field at a higher 
zoom factor. 
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Fig. 12.— Orbital solution for the 051(f) SB1 system MT734 based on 39 WIRO spectra 
from 2007-2011. 
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Fig. 13. — A time sequence of spectra for MT734, ordered by phase, in the spectral region 
around He I A5876 (left panel) and Ha (right panel). 



